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Difference Equation for Superrad lance* 
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We study the evolution of a completely excited system of N two-level 
atoms, distributed over a large region. Interacting with all modes of 
radiation field. We pay special attention to the distinction between r- 
conservlng (RC) and r-nonconservlng (RNC) processes. Considering the 
number of photons emitted as the discrete Independent variable, the evo- 
lution Is described by a partial difference equation. Numerical solution 
of this equation shows the transition from RNC dominance at the beginning 
to RC dominance later. This is also a transition from Incoherent to cohe- 
rent emlssi'~ of radiation. 
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In 1954 R. H. Dlcke^ discussed the spontaneous emission o£ radiation 

Iron jn excited system of N Identical two-level atoms. By considering the 

entire collection of atoms as a single quantum mechanical system, he found 

that under certain conditions the individual atoms cooperate to emit radia- 
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tion at a rate proportional to N which is much greater than their incohe- 
rent emission rate. This phenomenon is called superiradlance. The quantum 
number that plays the central role in the description of such system is the 
so-called cooperation number r which can be any integer between 0 and N/2, 
assuming N to be an even number. 
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In the vast literature on superradiance • most of the investigation 
have been limited to cases in which the emiting atoms either are confined 
to a region smaller than the wavelength of emitted radiation or are able to 
couple with only one radiation mode. In these situations r is conserved; 
which permits important mathematical simplifications. Current attention in 
this field is being focused on the more difficult, but more realistic, 
problem of N two-level atoms, distributed over a space of dimension greater 
than the radiation wavelength, interacting with all the modes of radiation 
field. In this situation r is no longer conserved. However, the effects 
of r-nonconservlng processes have not received sufficient attention in the 
literature. 
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Recently, we have developed a simple and systematic diagrammatic 
technique for calculating the matrix elements of collective operators; which 
gives the transition probabilities of both r-conserving (RC) and r-noncon- 
serving (RNC) processes as follows: 

T|^(r,m) « (r+m) (r-mfl) , (1) 
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( 2 ) 
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where k Is the wavevector of the photon emitted and Is a specific mode 
which may be the wavevector of the exciting pulse or the end-fire mode; m 
is 1/2 of the difference of the numbers of atoms in the excited state and 
the ground state, is Independent of r as long as r is very close 

to N/2; but it depends on the geometry of the sample. For a circular 
cylinder of length £ and radius a» it is known to be 






(3) 


where (k-kj^)^ p are the transversal and axial components of (k-kj^) and Jj^ 
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is the Bessel function of order 1. Rehler and Eberly have carried out 
the numerical calculation of which shows a very sharp peak along 

the direction of This implies that a photon emitted through RC process 
will be predominantly along the direction of while RNC process will 
give radiation along all other random directions. 

When (1) and (2) are Integrated over all possible directions, we obtain 

T^^(r,m) « (rfro) (r-iiri-l)p , (A) 


T^*^(r,m) “ [ (r+m) (r+m-l)/2rl (1 - w) , (5) 

where p has the Siime meaning as given by Ref. 6 and can have values much 
smaller than 1. 

Suppose we start with a completely excited system v.lth r ■» m ■ N/2. 
Then at the beginning photons will be emitted predominantly through RNC 
processes and, therefore, along random directions. For each photon emitted 
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through RNC process, both r and m decrease by 1; which does not Improve the 
chance of RC processes. On the other hand, for each photon emitted through 
RC process, r remains the same while ra decreases by 1; which will improve 
the chance of RC processes. Hence, the more photons emitted along 
direction, the higher the probability of additional photons emitted along 
this direction. The RC and RNC processes will be about equally important 
when (r-m)jj is of the order of 1. 

In this article, we are most interested in studying the transition 
from RNC dominance to RC dominance. To study the evolution of a system of N 
two-level atoms from a completely excited state, we will replace the time 
as the continues independent variable by the number of photons emitted as 
the discrete independent variable. Then our "dynamic equation" will be a 
difference equation. 

Let P(r,m) be the probability that our system has cooperation number 
r after N/2-m photons have been emitted. Then (4) and (5) gives the follow- 
ing partial difference equation for P(r,m) 


P(r,m) 






with the initial condition P(N/2,N/2) » 1 and the obvious restrictions: 

0 ^ T ^ N/2 and -r m ^ r. We have defined a = v/(l-p) which is essentially 
the same as s wher». }> 1$ 'i’ich less than 1. In obtaining (6) , we have ignored 
the possibility that a photon may be absorbed by the atomic system and the 
possibility that r may increase. The latter is justified in Ref. 5. 

It is convenient to define two new variables 


s H N/2 - r , t H r - m (7) 

and let Q(s,t) = P(r,m). Assuming that both s and t are much less than N, 
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(6) becomes 

ii 

Q(s,t) - [to/(l+to)]Q(s,t-l) + [1 + (t+l)o]"^(}(8-l,t) (8) 


with Initial condition Q(0»0) 1. 


By direct observation, the solution to (8) can be seen to be 


Q(s,t) 


* t+1 T r t T 

I n (l+io)"®i II (l+io)“^ tt 

_8 i»l JL^“1 


(9) 


where 6^"0, 1, 2, and the summation is over all possible choices of 

8^ such that 

t+1 

I 8. - s . (10) 

i-1 

The result of t'tls summation can be expressed as a contour integral as follows 

z = n (1 + iff) ^i 

8 i-1 

t+1 

z“®“^ n [1 + (l+ia)"^z + (l+ia)“V + ••• ] dz . (11) 
i-1 

A good approximation for Z can be obtained by integrating through a saddle 
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point on the real axis, the so-called steepest descent method , as follows: 



Z 




(l_:<^lo).. 
(1 - A + lo) 



t+1 

I 


i=l 


(1 + jq)A 
(1 - A + ia)2 _ 


( 12 ) 


where z - A is the location of the saddle point to be determined by the 
equation 


A ^ (1 - A + ia)“^ » s+1 
i=l 


( 13 ) 
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Substituting (12) Into (9) and using Stirling's formula for tl, ve obtain 


Q(s,t) « A”®(l+(t+l)o]| II (1-A+io) 

Li“l 

X e-' 


rt+1 


(n-t-)A 1 

JLl-1 (l-A+lo)^J 


- 1/2 


(14) 


Replacing s and t by the original variables r and we have a good 
approximation for P(r,m). Ue want to know what cooperation number has the 
highest probability after a certain number of photons have been emitted; 
i.e., we want to locate the maximum of P(r,m) with respect to r for fixed m» 
We will defined 


X = (N/2 - m)o , y = (N/2 - r)o , (15) 

consider x and y as continues variables, and replace all summations encoun*- 
tered by integrations. Then, from (13) and (14), we have 

A ln[l + (x-y)/(l-A)l - y , (16) 

lnP(y,x) “ ln(l+x-y) - (y/o+l/2)lnA + (l/2)lnt (x-y)/a| 

{ ■ 

+ [ (l-A+x-y)ln(l“A+x-y) + (l-A)ln(l-A) + (x-y)ln(x--y) ]/a 



Letting dlnP(y,x)/3y = 0 and using (16), we obtain 

l/(l+x-y) + l/[2(x-y)l + (l/a)ln[A(x-y)/(l-A+x-y)] 

+ (l-Hc-y)tA^(x-y) - (l-A)^(l+x-y)y]/2[A^(x-y) + (1-A) (1-A+x-y) 

- 0 . 


( 18 ) 
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Simultaneous numerical solution of (16) and (18) gives y as a function 
of X whlcht In turn, gives the most probable cooperation number as a function 
of the number of photons emitted. The results for five different values of 
a are presented In Fig. 1. A slop of 1 of these curves Indicates pure RNC 
processes and a slop of 0 Indicates pure RC processes. 

The author Is grateful to Mr. T. L. Kuo for his help with the numerical 


calculation. 
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FIGURE CAPTIONS 


Fig< 1* Drop In cooperation munber (N/2**r) as a function of the nunber of 
photons emitted (N/2-m) with a as a parameter. 
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